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Certain features associated with the symmetry reduction of the vacuum 
^. | Einstein equations by two commuting, space-like Killing vector fields are stud- 

ied. In particular, the discussion encompasses the equations for the Gowdy 
T 3 cosmology and cylindrical gravitational waves. We first point out a re- 
lation between the SL(2,R) (or SO (3)) a and principal chiral models, and 
then show that the reduced Einstein equations can be obtained from a di- 
mensional reduction of the standard SL(2, R) a-model in three dimensions. 
The reduced equations can also be derived from the action of a 'generalized' 
two dimensional SL(2,R) <r-model with a time dependent constraint. We 
give a Hamiltonian formulation of this action, and show that the Hamiltonian 
evolution equations for certain phase space variables are those of a certain 
generalization of the principal chiral model. Using these Hamiltonian equa- 
tions, we give a prescription for obtaining an infinite set of constants of motion 
explicitly as functionals of the space-time metric variables. 
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I. INTRODUCTION 

A central mathematical problem in non-linear field theories of physical interest is that of 
gaining 'control' over the space of solutions. One strategy is to map the theory of interest 
in to 'standard' models which have been well-analyzed. Another is to find the 'full set' of 
conserved quantities associated with field equations. Both these strategies have proved to 
be useful in two-dimensional field theories, where a number of integrable models have been 
known to exist for some time [1|,Q] . Now, in presence of two commuting, space-like Killing 
vector fields, the four dimensional, vacuum Einstein theory reduces to a two dimensional 
field theory. Hence, it is natural to explore the structure of this sector of general relativity 
using techniques that have been successful in standard integrable models. 

If the Killing fields under consideration are hypersurface orthogonal, the theory has only 
one local degree of freedom, and the field equations can in effect be reduced to a rather 
simple linear partial differential equation. The model is then exactly soluble and one can 
explicitly write down the full set of conserved quantities 0-§| . On the other hand, when the 
Killing fields fail to be hypersurface orthogonal, the theory is both more difficult and more 
interesting: Now, there are two local degrees of freedom which remain coupled non-linearly. 
The space-times in this category include the Gowdy cosmological models || and cylindrical 
gravitational waves (see |7|] and references therein). In this paper, we shall focus on these 
models. 

Various aspects of the underlying equations have been studied in detail over the years 
(see, e.g., |f- |Ll| ). The results have led to a wealth of insights in to the structure of the 
space of solutions of the reduced theory which, in turn, have made these 'midi-superspaces' 
convenient tools to probe physical issues that arise in the full theory, such as quantization 
and cosmic censorship. Quantization of the one polarization cylindrical waves was studied 
in |7j and more recently, in |IS| . The second analysis was subsequently used to show that, 
in this model, there exist surprisingly large quantum gravity effects which undermine the 
validity of the classical and semi-classical theories in unexpected ways W7\ . More recently, 



a quantization of the more interesting two-polarization gravitational waves has also been 
achieved |TS| and its physical implications are now being explored. Finally, the models 



have proved to be useful in analyzing cosmic censorship in cosmological [T^] as well as 
asymptotically flat contexts [|2D| . 

The investigations of |8|-|15|j have also shed light on the mathematical structure of the 
reduced Einstein equations. In particular, the presence of the infinite dimensional Geroch 
group which acts transitively on the space of solutions indicates that the system is in some 
sense exactly integrable. The application of the inverse scattering methods has also led to 
results supporting this view point. More recently there has been work on obtaining constants 



of motion for the reduced Einstein equations with open spatial topology in |21] and closed 
spatial topology in |plfl . However, it is fair to say that we do not yet know if there exists 
a complete set of constants of motion; unlike in the standard, two-dimensional integrable 
models MM, an analysis involving, e.g., bi-symplectic structures has not been carried out. 
Hence it is of interest to seek new avenues to this problem. 

The purpose of this paper is two fold: i) to point out new relations between the re- 
duced Einstein system and the other well-known models which might be useful in future 
investigations; and, ii) to propose a new strategy to address the problem of finding explicit 



constants of motion. More precisely, we will study the reduced Einstein equations in the 
metric variables, further develop the interplay between these equations and the standard 
and generalized cr-models in three and two dimensions respectively, and, using this relation, 
provide a new avenue for obtaining an infinite set of constants of motion. 

The plan of the paper is as follows. In Sec. II we recall the reduced Einstein equations 
for the Gowdy T 3 cosmology and cylindrical gravitational waves. In Sec. Ill, we first 
describe a relation between the SL(2, R) or 50(3) cr-models and the principal chiral model 
in any space-time dimension, and then show that the reduced Einstein system is a sub- 
system of the standard, three dimensional SL(2, R) cr-model. In Sec. IV we elaborate on 
the relation between a generalized cr-model in two dimensions, and present a Hamiltonian 
formulation and associated evolution equations. This framework is utilized in Sec. V to give 
a prescription for obtaining an infinite set of constants of motion. 

II. PRELIMINARIES: REDUCED EINSTEIN EQUATIONS 

Recall that solutions to the vacuum Einstein equations with two commuting, space-like, 
Killing vectors fall into classes of well known space-times. These include plane, cylindrical |7j 
and 'toroidal' [E3| gravitational waves for non-compact spatial topology and the T 3 , S 1 x S 2 



and S 3 Gowdy cosmological models for compact spatial topology |J. For definiteness, we 
will consider one representative from each class, namely the T 3 cosmological model and 
cylindrical gravitational waves. Although this material is well-known, it is included here 
because the final results will serve as points of departure for our main analysis. 

A. Two polarization Gowdy model 

The T 3 Gowdy model has the space-time metric 

ds 2 = e 2A {-dt 2 + d8 2 ) + g ab dx a dx b , (2.1) 

where A = A(t, 8), x a = (x 1 , x 2 ), and the 2x2 metric 

„ /coshW^ + cos$ sinhW 7 sin$ sinhiy \ , > 

~~ V sin$ sinhW 7 coshiy - cos$ sinhW 7 J ^ ' ' 



is parametrized by the three functions R(t, 8), $(£, 8) and W(t, 8). The metric ( [2.1|) has two 
commuting space-like Killing vector fields 

/ ° ' ^ (£) , (2.3) 



y dx 1 J \ dx 2 

with the additional restriction that each of their orbits has the topology of a two torus. 
The t, 8 coordinates have ranges < t < oo and < 8 < 2tt. With these conventions, the 
metric is that of the 3-torus Gowdy cosmology ||. The vacuum Einstein equations in the 
Gowdy gauge, defined by setting R(t, 8) = t, give the two coupled two-dimensional evolution 
equations 



W + -W-W" + sinW cosW ($ /2 - $ 2 ) = , (2.4) 

L 

$ + 1$ _ $" + 2 ^^ ($iy - <W) = , (2.5) 

£ sinhW 7 

for W(t, 9) and $(£,#), and the two 'constraint' equations 

A + j - - [ W 2 + W 2 + smh 2 W ($ 2 + $' 2 ) ] = , (2.6) 

A' - - (WW + smh 2 W $$') = . (2.7) 

Note that the evolution equations involve only W(t, 9) and <&(t,9) and, given a solution 
to these equations, the field A(t, 9) can be obtained by a simple integration of the con- 
straint equations. (Eqns. ( |2.4| ), ( |2.5| ) also serve as the consistency conditions for the two 
constraints.) 

However, because the spatial topology is compact, there is a subtlety: We can solve for 
A consistently only if the integral over the circle of the expression defining A' in Eqn. ( |2.7|) 
vanishes. [] This leads to a global constraint relating the two fields W(t, 9) and $(£,#), 
namely, 



r d9 (WW + sinh 2 H/ $$') = , (2.8) 

Jo 



f-271 

/O 

which is preserved by the evolution equations. Thus, to obtain a solution to the Gowdy 



model, we need to solve ( |2.8| ) and the evolution equations Q2.4J ) and (|2.5|). However, as we 
will see, in the Hamiltonian framework, the global constraint generates just the rigid rotation 
around the circle. Therefore, it has rather simple Poisson brackets with other constraints 
and constants of motion. Hence, the heart of the problem lies in the evolution equations. 
We shall initially focus on these, and incorporate the global constraint at the end. 
The equations, ( |2.4| - [2.5|) , may be derived from the action 



S G (W, $) = - f dtd9 t[W 2 -W 2 + sinh 2 H/($ 2 - $' 2 ) ] 

= I J dtd9 t v^ V ab G AB (Y) d a Y A d b Y B , (2.9) 

where Y 1 = W, Y 2 = $, r] ab = diag(-,+), a, 6, ■ ■ ■ = t,9, and G AB (Y) dY A dY B = 
dW 2 + sinh 2 ^ <i$ 2 is the unit hyperboloid metric. The Gowdy model with one polarization 
is obtained by setting $ = in the above equations. This leads to a linear evolution equation 
for W, and hence exact solvability. 

B. Cylindrical gravitational waves 

The space-time metric for cylindrical gravitational waves can be put into the general 
form (|2.1| - |2.2|) above. The line element is 



l We thank the referee for drawing our attention to this point. 



ds 2 = e 2A {-dt 2 + dr 2 ) + g ab dx a dx b . (2.10) 

The spatial topology is now R 3 , the orbits of the two Killing vector fields has topology 
S 1 x R, and all the metric functions depend only on t and r. Now the relevant gauge fixing 



is R(t,r) = r. The vacuum Einstein equations are very similar to (p74" -275); all the t and 
r factors are interchanged with appropriate sign changes. The evolution and constraint 
equations are 

W--W-W + sinhiy coshiy ($' 2 - $ 2 ) = , (2.11) 

$ _ i $' _ $" + 2 ($W - &W) = , (2.12) 

r sinhVT 

and 

A' + -!- - J [ W^ 2 + I^' 2 + sinh 2 I^ ($ 2 + $' 2 ) ] = , (2.13) 

i - - (WW + smh 2 W $$') = , (2.14) 

where now ' = d/dr. Since the spatial topology is now R 1 , there is no longer a global 
constraint. We can therefore focus simply on the evolution equations. These follow from 
the action 

S C (W, $) = - fdtdr r[W 2 -W' 2 + sinh 2 W($ 2 - $ /2 ) ] . (2.15) 



The only difference in the Lagrangian densities of (2.9) and (|2.15|) is the interchange of the 
t and r factors. 

The awkward feature of these actions is the explicit factor of t or r in the respective 
integrands. Without these factors, the actions would be those of the usual two-dimensional 
non-linear cr-model, and standard methods for obtaining complete sets of conserved currents 
would be applicable. In the present case, these procedures have to be modified. 

We will discuss two possible avenues. The first is to embed the two dimensional model 
( [2.9|) in a standard three dimensional, flat space non-linear cr-model, which is free of the t (or 
r) factor. (The two-dimensional model, with its t (or r) factor is recovered as a dimensional 
reduction of the standard three dimensional model.) In this approach, the solutions of the 
reduced Einstein theory constitute a specific (symmetry-reduced) sector of the solution space 
of the three dimensional cr-model. The second avenue is to remain in two dimensions and 
handle the t (or r) factor by imposing a t (or r) dependent constraint on a free theory of 
three scalar fields. We will explore these avenues in sections III and IV respectively. 

III. RELATION TO THREE DIMENSIONAL a AND CHIRAL MODELS 

In Section III. A, we point out a relation between SL(2, R) or SO (3) cr-models and prin- 
cipal chiral models for these groups which exists in any space-time dimension but which 
appears not to have been noted in the literature. In Section III.B, we return to the reduced 
Einstein system and show that its space of solutions can be naturally embedded in to the 



space of solutions of the SL(2, R) cr-model in three space-time dimensions. Thus, the three 
models are closely related; the reduced Einstein theory is a sub-case of the other two mod- 
els. In particular, any method of showing integrability or extracting conserved quantities in 
three dimensional a or chiral models is directly applicable to the reduced Einstein system. 

A. <7 model and flat connections 

For concreteness let us consider the SL(2, R) non-linear cr-model in n-dimensional 
Minkowski space-time; with obvious sign changes, our discussion goes through step by step 
in the SO (3) case and for Euclidean space-time signature. The model has the action 

S n (X, A) = -i J oTx V=V [ r, ab daX l d h Xig l3 + A (gijX'X* + 1) ] , (3.1) 

where r] ab is the flat Lorentzian metric of signature (— , +,•■■,+), g^ is the Cartan metric 
for SL(2, R) of signature diag(+ H — ), and X(x) is a Lagrange multiplier. The equations of 
motion are: 

□X ! = rfd a d h X i = XX 1 (3.2) 

gijX { X j + 1 = 0. (3.3) 

To see the relation to the principal chiral model, define a Lie algebra valued 1-form by 

A\ = -2f jk X*d a X\ (3.4) 

where /*- A are the SL(2, R) structure constants. Then it follows from the equation of motion 
(U) that 

d a A\ = . (3.5) 

Furthermore, using X l d a Xi = 0, which follows from the constraint (|3.3|) , it is easy to show 
that 

J& = d a A\ - d h A\ + f jk A{A k b = . (3.6) 

Thus, the cr-model equations imply ( |3.5|) and (|3.6| ), which are the equations of the principal 
chiral model. The last step makes a crucial use of the three dimensionality of our Lie groups 
(i.e., the specific form of the structure constants of SL(2, R) (or 5*0(3))). Hence, while this 
relation holds in any space-time dimension, it does not hold for general groups. 

Let us now ask if this map from the the space of solutions of the non-linear sigma model 
to that of the principal chiral model is surjective. Can we derive ( p.2| - ^.3D from the chiral 
model equations? First, F ab = implies that A a has the form A a = g~ x d a g for g e SL(2, R). 
Now, a standard parametrization of g is 

( a + X 2 X x + X 3 \ 
9={ x i_ x s a _ X 2 ) , (3.7) 

with the unit determinant condition a 2 — gijX l X^ = 1. This gives 



A a = g-'dag = 2 (ad a X - Xd a a) - 2 [X, d a X] , (3.8) 

where X = X l Ti is an element of sl(2, R) with generators r». Secondly, d a A a = leads to 
the equation 

[X ) UX] = aUX -XUa . (3.9) 



The unit determinant condition implies that the cr-model constraint ( [3.3|) is satisfied if and 
only if a = 0. In this case, the commutator [X, OX] vanishes which in turn implies that 
OX 1 = XX 1 for some function X(x). Thus, the cr-model equations hold if and only if the 
SL(2, R) matrices g are trace free. To summarize, the space of solutions to the cr-model 
is naturally embedded in the space of solutions to the principal chiral model. The image 
consists of flat connections A a given by A a = g~ 1 d a g where g is trace-free (i.e., where a = 
in the parametrization ( |3.7| )). We will return to this result in the next subsection. 

To conclude this discussion, we note that there is also a Hamiltonian version of the 
derivation of the SL(2, R) principal chiral model from the corresponding cr-model, which 
is presented below for a more general time dependent constraint X l Xi + 1 = 0. As we will 
see, this is the type of constraint relevant for the reduced Einstein equations if one wants to 
work with two dimensional models. 



B. Reduced theory from three dimensions 

Let us now specialize to three-dimensional space-times with topology S 1 x R 2 , equipped 
with a flat metric r\ abn given by 

r] ab dx a dx b = -dr 2 + dx 2 + d6 2 . (3.10) 

Consider the action 

S' 3 [Y] = ~\ I drdxdO y^j v ab G AB (Y) d a Y A d b Y B , (3.11) 

I JM 

where, as before A = 1,2 but now Y A = Y a (t,x,9). As is well-known, this action yields 
the same equations of motion as Ss(X, A) of ( |3.1|) . Note, however, that, in contrast with the 
reduced Einstein action ( |2.9| ), there is no factor of r in the integrand of ( |3.11| ). However, 



we will now show that (|2.9| ) does result from (|3.11|) by a symmetry reduction. Thus, the 
reduced Einstein model is contained in the three dimensional, standard, non-linear cr-model. 
To perform the required reduction to two dimensions, we proceed in two steps. First, let 
us make the change of coordinates 

t = t costu/, x = t sinhy (3-12) 

in the r — x plane, which casts the metric in the form 

rj ab dx a dx b = -dt 2 + t 2 dy 2 + d9 2 . (3.13) 

Second, let us require that the Lie derivative of the field variables Y A (t,y,9) along the 
boost-Killing vector field (d/dy) a of r] ab vanish. (In the chosen coordinates this means 
Y A = Y A (t,6).) Then, we have: 



S' 3 [Y] — S' 2 [Y] = -~ JdyJdtdO t V ab G AB (Y) d a Y A d b Y B = S G [Y] (/ dy^j , (3.14) 



which is the Einstein action ( |2.9| ) for Gowdy models, multiplied by the constant / dy. 

Thus, we arrive at the following conclusions. As we have just shown, the space of solutions 
Sq to the reduced Einstein equations (i.e. to the Gowdy model) is naturally embedded in 



the space S a of solutions to the cr-model resulting from the action S' 3 [Y] of (|3.11| ), or S3 [A, A] 
of ( |3.1|) . Second, we saw in section III. A that the space S a is in turn embedded into the 
space of solutions S c h to the principal chiral model. The image of the first map consists 
of fields Y A (t, 0) which are ^/-independent, and satisfy the global constraint (2.8), while 
the image of the second map consists of A a = g~ 1 d a g for which the g of ( |3.7| ) is trace- 
free. Therefore, any method for finding conserved quantities for the standard non-linear 
cr-model or the principal chiral model in three dimensions yields, via restriction, a method 
for finding conserved quantities for the reduced Einstein equations. However, in practice, 
this simplification is not directly useful since the standard techniques for finding conserved 
currents for a and chiral models are tailored to two space-time dimensions. Nonetheless, the 
existence of embeddings is conceptually interesting and may have a more general application. 
For example, results on integrability, asymptotic forms of solutions, and quantization of 
the standard three dimensional a and chiral models will also apply to Gowdy models and 
cylindrical waves by restriction. Indeed, this is essentially the procedure used to (find the 
complete set of constants of motion and) quantize one polarization cylindrical |7|JT6| , toroidal 
|5|] waves and Gowdy models [TJ|]. 



IV. RELATION TO TWO DIMENSIONAL GENERALIZED a AND CHIRAL 

MODELS 

Let us now explore the relation between the reduced Einstein system and two dimensional 
models. In three space-time dimensions, we could restrict ourselves to the standard a 
and chiral models. In two dimensions, on the other hand, we will have to allow certain 
generalizations. 

More precisely, in Sec. IV. A, we will show that the reduced Einstein equations can be 
derived from a two dimensional generalized cr-model action, where, however, the constraint 
is time (or space) dependent. In Sec. IV. B, we will perform a Legendre transform and 
show that the Hamiltonian evolution equations for certain variables give a generalized chiral 
model. This latter result is a Hamiltonian version of that given in Sec. III. A above, suitably 
generalized to incorporate the space or time dependent cr-model constraint. 

A. Time dependent <r model 

Let the space-time topology be S 1 x R and let r] a b be the obvious flat metric with 
signature (—,+). Consider the following two-dimensional action for three free scalar fields 

X*(i = 1,2,3): 

S[X, M = ~\ J dtd6 [y/^j r? h g l3 d a X l d b X ] - A (g tJ X l X j + t)] , (4.1) 



where g^ =diag(+ H — ). The variation of S with respect to A gives the time dependent 
constraint 

d := g ij X i X j + t = X ■ X + t = . (4.2) 

This constraint is solved by setting 

X 1 = Vt cos$ sinhW, X 2 = Vi sin$ sinhW, X 3 = Vt coshW , (4.3) 

where the fields W(t, 9) and $(£, 6) are unrestricted. These conditions, when substituted in 
to the action (|4.1|) , give the reduced action 

S R [W, $\ = \l dtd9 I t [ W 2 - W' 2 + sinh 2 iy ($ 2 - $' 2 ) ] - — | . (4.4) 



The integrand here differs from that of Sg of fl2.9|) only by an additive field independent 
term, and therefore leads to the same evolution equations as ( |2.9| ). It is also straightforward 
to verify directly that, as is expected on general grounds, the equations of motion following 



from S[X, A], together with ( |OD , lead to the reduced Einstein equations (|2.4|- |2.5|) . 



B. Hamiltonian formulation and generalized chiral model 



We will now show that the Legendre transform of the action (|4.1| ) leads to "chiral model- 
like" equations, where however the curvature is not flat. The phase space form of the action 

is: 

S[X, X, P] = J dtdO I P t X l - i [g^P^ + g^X[X ! 3 - A {g^X l X j + t) ] } . (4.5) 

where Pi are the momenta conjugate to X % . Following now the Dirac prescription, we require 
that the constraint C\ of (|4.2j) be preserved in time: 



f/r ' {Cu H} + ^ = 2X-P + l = 0, (4.6) 



dt ' ' dt 

where 

1 fl-K 1 /■27T 



h = - do ^{PiPj + x[x' 3 ) + - / de \d 

2 Jo J 2 Jo 

= H + l Fffl\C x (4.7) 

2 Jo 

is the Hamiltonian identified from ( |4.5| ). This gives the secondary constraint 

C 2 :=2X-P+1 = 0. (4.8) 

The pair of constraints C\ = and C<i = are of second class. The constraint Ci is solved 
by setting 



P 1 = — -= cos<3> sinhW 7 - \ft sin$ sinhW 7 $ + \ft cos$ coshW 7 W 
2\Jt 

P 2 = — f sin$ sinhW 7 + \ft cos$ sinhVT $ + \Tt sin$ coshW/ W 
2y/i 

P 3 = F coshW/ - Vt sinhW W 

2y/i 



(4.9) 



and using (|4.3| ). As we did for C\, let us require that the secondary constraint C 2 be preserved 
in time. This will lead either to a new constraint, or to a condition on the lagrange multiplier 
A. It turns out that the latter is the case : 



<K ' 2 {C 2 ,H} + ^ = {C 2 ,H } + 4tX=0 



dt 



dt 



which gives A = — {C 2 , H Q }/4t and 

H = H -± : T de{C 2 ,H Q }C x . 

at Jo 



(4.10) 



(4.11) 



This is the (first class) Hamiltonian we must use to derive evolution equations. The above 
procedure is equivalent to using the Dirac brackets for the unmodified Hamiltonian Hq. To 
summarize, a Hamiltonian description of the evolution equations of the Gowdy model can 
be given as follows. The phase space consists of canonical pairs (X\Pi). There are two 
second class constraints and the Hamiltonian is given by Q4.ll ). If we initially satisfy the 



two constraints, they continue to be satisfied by the Hamiltonian flow. 

Recall however that the equations governing the Gowdy model also includes a global 
constraint ( |2.£| ). It is straightforward to verify that this constraint generates the rigid 
diffeomorphism on S 1 . Hence it weakly Poisson-commutes with C\ and C 2 and strongly 
Poisson-commutes with the Hamiltonian. Therefore, it is a single (global) first class con- 
straint within the Hamiltonian description arrived at above and we can just carry it along 
in the analysis. Constants of motion of the system without the global constraint will provide 
physical constants of motion simply by restriction to the part of the phase space where the 
global constraint is satisfied provided they are invariant under rigid diffeorphisms of S 1 . 

To go to a generalized chiral model, it is convenient to replace (X 1 , Pj) pairs by new 
variables 



L i = f i jk X j P k and J i = f i ik X i X 



k! 



)>-' 



(4.12) 



where f l j k are the SL(2, R) structure constants. We will now argue that there is no loss of 
information in using L l , J 1 in place of X 1 , P % satisfying the two constraints (|4.2j ) and ( [1.6| ). 
First, it is easy to verify that both L l and J 1 are space-like in the internal directions (but 
can vanish). Hence, J x L is time-like. When this last vector is non-zero, one can easily 
recover X % and P l from L l and J 1 explicitly: 



X 



Vi 



J xL 



— J x L , 



P = ± (-J x L 4 (J x L)x L 

Vi JxL V2 



(4.13) 



where 



10 



J x L 



-(J x L) • (j x L) . (4.14) 



The choice of the sign in the two equations can be fixed by first requiring that X % be future- 
directed and then using ( |4.6| ). Since J x L is time- like, its norm can vanish only when the 
vector itself vanishes. In this case, we can not simply use explicit formulas ( |4.13| ). However, 
a more detailed, systematic analysis shows that X 1 , P % can be again recovered from L l , J 1 . 
L % and J 1 satisfy simple Poisson bracket relations 

{L i (t,9),L j (t,e'} = f ij k L k (t,9)6(9,e'), 

{J t (t,9),J 3 (t,9')} = 0, 

{Li{t,6), Jfaff)} = 4- fc J fc M) S(9,9')+ gij 8'(9,9') . (4.15) 

Furthermore, 

{L%9),C x {t,9')} = {U(t,9),C 2 (t,9')} = (4.16) 

{ J% 9), C t (t, 9')} = 0, { J l (t, 9), C 2 {t, 9')} = U% 9) 5(9, 9') . (4.17) 

(The indices i,j,--- are lowered (raised) using the SL(2,R) metric g^.) The evolution 
equations of the SL(2, R) variables (L l , J 1 ) are: 

U = {L\H} = J*', (4.18) 

where, as before, the prime denotes the derivative with respect to 9, and 

X = { J\ H } - {J\ d} Q) {C 2 , H} = {J\ H } 

= V' - if jk X j 'P k . (4.19) 

The last equation can be further simplified by noting that 

pi ^ h xi + \ fl > kX ' Lk > (4 - 20) 

« -7T + - PikL j J k , (4.21) 

2t t 3k K ' 

where ~ denotes equality modulo the constraints C\ and C 2 - This finally gives 



\W + \ 



ji = u< _ -f jk L j J k + -J\ (4.22) 



We are now ready to put these Hamiltonian equations in to a chiral model-like form. 
Define the matrices 

A := 2V T i , A x := 2JV; , (4.23) 

where 

11 



1 /O 1 



Tl = 2Vl 



1 (\ 



T2 = 2V0 



T3 = 2 



(4.24) 



_f k 

Jij T k ■ 



are the generators of SL(2, R) satisfying the relations 

1 

[^ T j] = fij kr k , TiT 3 

Multiplying the evolution equations Q4.18 ) and ( |4.22| ) by r* leads to 

d A - dxA x = 
d Q A 1 - d x A + [A , A ± ] = (l - 1) [A , A,] + ^ , 



(4.25) 

(4.26) 
(4.27) 



where 0, 1 = t, 9; t > and < 9 < 2n. This is the desired chiral model-like form of the 
Einstein evolution equations ( |2.4| - |2.5| ). 

So far, we have considered the Gowdy model equations where the relevant a-model 
constraint is time dependent. The analysis of cylindrical waves is completely analogous. 
The evolution equations (|2.11| - fET2D may be derived from an action like (|4.1| ), but with r 
replacing t in the constraint. After performing a Hamiltonian decomposition, and following 
steps similar to the above, we find that the cylindrical wave evolution equations for the 
variables L, and J; lead to the chiral model-like equations 



d A - d x A x 
d Q A x - dxA + [A Q , Ai] 







[^o, Ax] + 



A, 



(4.28) 
(4.29) 



Note that these are different from the Gowdy model equations fl4.26 - 4.27 ) only in the replace- 
ment 1/t, Ai/t — > 1/r, Ao/r. In addition, now, there is no analog of the global constraint 
( p.8|) to worry about. 



V. CONSTANTS OF MOTION 

For the chiral model in two space-time dimensions, there exist standard techniques to 
construct an infinite hierarchy of conserved currents. In this section, we will show that one 
of them admits a suitable extension which is applicable to the generalized chiral model of 
Sec. IV. 



A. Gowdy cosmology 



The idea is to use the Hamiltonian equations ( f4.26 - 4~2^| ) to extract constants of motion 
for the reduced Einstein equations. The method is a modification of one given by Brezin et. 
al. (BIZZ) for obtaining constants of motion of the principal chiral model p5| . Let us first 
rewrite the evolution equation (|4.27 ) in a convenient form: 



01 



Ax + il-^daA, 



(5.1) 
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In the standard chiral model, the right side vanishes. The divergence-free condition (|4.26|) 
on A a is the same in the two cases. So, modification of the (BIZZ) procedure of |25j is 
necessary to handle the right side of (|5.1|). 

As for the principal chiral model, (|4.26|) is already a conservation equation, so we can 
write the first conserved current as 



J® = DaX {0) 



(5.2) 



where D a = d a + A a and A^ -* is any constant matrix. Denote the corresponding conserved 
charge by Q : 



Q 



(i) 



2- 



de j, 



(i) 



2?r 



de 



An, A<°> 



(5.3) 



Since J^ — (Q^ /2n)d a t =: J^> — d a (p^ is a smooth conserved current with zero charge, it 
is dual of the gradient of a smooth field. Thus, there exists a smooth, matrix-valued field 
A* 1 ) such that 



This equation determines A^) up to an additive constant matrix: 

a«(M = - f e dff 4 1 \e',t) + ^-~[ t dt' jP{o,i!) . 
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(5.4) 



(5.5) 



Note that, in spite of the explicit appearance of 9 in the second term, the 'potential' A^ is 
smooth everywhere, including the point 6 = 0. It depends on A a , i.e., on the solution to the 
Gowdy model under consideration and the idea is to construct another conserved current 
from \^\ 

As in the (BIZZ) procedure of |2"5| , let us introduce a fiducial current, 



< 2) : = D a \V , 
and compute its divergence 

V ab d a Ki 2) = ri ab D b d a \ {1) = e ab D a (j fe (1) - d b 4>^) 
= e ab D a D b \^-e ab \A a ,d b ^' 



(5.6) 



>oi,A<°> 



A a ,d b ^ 



Al ,M+(l- t )e ab \d a A b ,\W 



— e 



ab 



A a ,d b ^ 



e ah d n 



(1 - t)A b , A^l + e ab S° a \A b , A<°>1 + U, A<°> 



ab 



A a ,d b ^} , (5.7) 



where the equation of motion ( |5.1| ) is used in the fourth equality. Collecting terms, we obtain 



V ab daK 



(2) 



e ab d n 



1 - t)A b , A (0) l + 2e ab 5° a \A b , A (0) l + e ab 5° a 



A h 



2n 



(5.8) 



Now the key observation is that the last two terms on the r.h.s. may be rewritten as 
divergences. Indeed, using (|5.4|), we obtain 
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e ab 5 Q „ 



A b , A<°>] = e ab 5° a (e b c d c X {1) + <9 6 0«) = V ab d a (X^5° b + 0«^) . 



(5.9) 



Similarly, [A a , Q^ 1 ' /2tt] is a conserved current. Using it in place of J^ = [A a , A^] in the 
steps that led us to Eqn. (|5.5|), we obtain 



afcrf) 



e ab 5 



A h 



2vr 



e ab 6° a (e 6 c 5 c A« + d b ^) = rj"»d a (X^6° b + ^8\ 



(5.10) 



Here A 1 - 1 - 1 is defined similarly to A*- 1 - 1 of Eq. ( |5.5j ) by 



\ (1 \9,t) 



d0' A (9',t),Q {1 >/2ir 



Q« 



2ty 



dt' 



U, 



Ai(0,f 



, Q (1) 



2vr 



(5.11) 



with 



A a ,QW/2n] replacing jW = U a , A<°> 



and 0W = Q^H/2n with 



Q« = / dd 
Jo 



2 IT 



A (9,t),Q^/2n 



Thus, we have shown that the current 

Jf ee K& - e a \l - t) [A, AW] - S° a (2A« + A«) - # (20« + 0«) 
is conserved. The corresponding conserved charge is 



(5.12) 



(5.13) 



Q 



(2) 
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d# j,S 2) ^,t) 



2?r 



dfl { D A (1) - (2A (1) + A (1) ) + (1 - t) [A u A (0) ] } , (5.14) 



with A^ and A^ are defined as in (|5.5|) and ( |5.11 ). 



Notice that Q^ = Q\ Vj and Q^ = Q\ V* are independent charges because the latter 
depends explicitly on t, Aq and A\, whereas the former depends only on Aq. These are in 
fact six different charges because we can "peel off" the matrices Tj. By inspection, these 
charges are invariant under rigid diffeomorphisms of S 1 . Hence, they Poisson-commute with 
the global constraint ( |2.8| ) and are physical constants of motion for the Gowdy model. 

The key question now is whether we can build a whole tower of conserved charges. 
A natural strategy is to use Poisson brackets between the Q\ and Q\ . The brackets 
between Q\ among themselves form the sl(2, R) Lie algebra; they do not yield new charges. 
Moreover, their brackets with Q^ just rotate the Q\ among themselves. On the other 
hand, the Poisson brackets {Q\ ,Q) } provides charges whose expressions have a higher 

(2) 

degree of non-locality than the Q\ themselves; the expressions contain one more nested 
integral. It is therefore reasonable to suppose that, generically, the Poisson bracket is a new 
conserved charge. It is straightforward to see that the continuation of this process gives 
conserved charges of successively higher degrees of non-locality, each of which appear to 
be functionally independent of the previous ones. Thus, the procedure yields an infinite 
tower of conserved quantities. Although we do not have a conclusive proof, the structure 
of non- locality suggests that all these charges are independent. Finally, by construction, all 
these charges Poisson-commute with the global constraint ( |2.8| ) and are therefore physical 
charges for the Gowdy model. 
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Let us summarize. The equations of motion following from the action (|4.1|) are the 
reduced Einstein equations. We performed a 1 + 1 -decomposition of this action to obtain 
the Hamiltonian theory, and then made a change of variables from (X l ,Pi) to A l a . The 
equations for A\ resemble those of the chiral model. Hence, it is possible to extend the 
standard procedure |25| to obtain two sets of conserved charges Q\ and Q\ • . A tower of 
new conserved quantities can be obtained by taking Poisson brackets. Finally, note that the 
constants of motion can be rewritten in terms of the original metric variables (W, $) and 
their space and time derivatives by using the solutions ( |4.3| ) and (|4.9|) of the constraints C\ 
and Ci- 



B. Cylindrical waves 

As noted before, the evolution equations for cylindrical waves and the T 3 Gowdy model 
are very similar. Therefore it is not surprising that the above approach to finding con- 
stants of motion applies to cylindrical waves. There are two main differences. First, the 
spatial topology is now R 3 rather than T 3 , whence the spatial sections of the reduced two- 
dimensional model are half-lines < r rather than circles. The triviality of this topology 
will simplify the analysis. The second difference is that we now need to impose boundary 
conditions on our dynamical variables. For simplicity, we will assume that, on an arbitrarily 
chosen initial Cauchy surface, the fields W and $ are of compact support. This assumption 
will ensure the convergence of various integrals but can be weakened in an obvious manner. 

As in the Gowdy model, the existence of the first conserved charge follows immediately 
from (|4.28| ). The conserved current has the same form as before: jfp = D a fi^°\ where fx^ is 
any constant, fiducial matrix. However, because the spatial topology of the reduced model 
is now trivial, the potential jjP~> is now defined simply by 



4 1] = Ca%» {1) , 

which determines //^ up to an additive constant matrix, 

(J.W = - f dr 1 4 1] (r',t) - f dt' J! (1) (0,O . 
jo Jo 

The second, fiducial current has the same form as J^ 

4 2) := D a ^ . 



(5.15) 



(5.16) 



(5.17) 



V ab da^ 



e ab D n J, 



Taking its divergence, we find 

e ab D a D b ^ 

--(l-r)e ab [d a A b ,^] + [A ,^ 

e ah d a [{l-r)A b ,n m 
V ab d a \(l-r)A c ,^°y 



Foi,fi> 



(i) 

a°b 
(0)" 



(5.18) 



This allows us to identify the second conserved current 
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4 2) = 4 2) - (1 - rK h [A b , „<°>] . (5.19) 

Thus, the triviality of spatial topology simplifies the analysis. The corresponding conserved 
charge is 

/*oo P oo 

Q {2) = dr 4 2) = dr {D of i^ + (1 - r) [A, //>)] } . (5.20) 

As before, further conserved quantities can be generated by taking Poisson brackets. 
Finally, in this case there is no global constraint to take care of. 



VI. DISCUSSION 

Let us summarize the main results. In Sec. Ill, we first showed that the standard, 
SL(2,R) non-linear a-model is embedded in the principal chiral model in any space-time 
dimension. We then showed that the symmetry reduced Einstein system is embedded in 
the a model in three space-time dimensions. Thus, there is a hierarchy and results from the 
three-dimensional a and chiral models can be taken over to the reduced Einstein system. 
This strategy is successful in the one-polarization case, i.e., the case when the two Killing 
fields are hypersurface-orthogonal. 

In Sec IV we worked in two space-time dimensions and recast the reduced Einstein 
model as a 'time dependent' <r-model, or, alternatively, 'generalized' chiral model. For 
the standard chiral model, there exist the so-called BIZZ procedure which enables one to 
construct a hierarchy of conserved charges. In Sec. V, we showed that the procedure can 
be appropriately modified to the generalized chiral model of Sec. V to obtain the analogs of 
the first two BIZZ conserved currents. Modifications were required because of two reasons: 
i) the two-dimensional space-time topology is S 1 x R, rather than R 2 ; and ii) the curvature 
of the generalized model does not vanish, but has a specific form. (To our knowledge, the 
complications that arise due to non-trivial topology have not been discussed in the literature, 
even in the standard BIZZ procedure.) Since the currents take values in the Lie-algebra of 
SX(2,i?), we obtain six conserved quantities, Q\ and Q\ , where the index i refers to the 
Lie-algebra. One can restrict this procedure to the well-understood, one-polarization case. 

(2) 

Even in this case, the second, non-local set of charges Q\ appears to be new, i.e., appears 
to be a non-trivial combination of the known charges. 

(2) 

Returning to the general case, Poisson brackets between the Q\ lead to new conserved 
charges and one can continue the process by taking Poisson brackets between the available 
charges. The functional form of these charges exhibit increasing non-locality, each step 
leading to an additional 'nested integral'. The initial- value problem for these equations is 
well-posed, and each charge appears to probe a different aspect of the functional form of 
the initial data. Therefore, it appears that the charges are all functionally independent. 
However, we do not have a definitive proof of this independence. For example, after a 
certain stage, the Poisson bracket may just yield a c-number, i.e., a constant functional on 
the space of solutions. In this case, the procedure would terminate. It may also happen that 
at a certain stage the conserved charges cease to be functionally differentiable, making it 
impossible to compute further Poisson brackets. However, the non-local structure of charges 
is reminiscent of the structure of the commutators in the Geroch group. Is there a close 
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relation between the two? If so, one would have an elegant avenue to show that the charges 
are indeed independent. 

Finally, we should emphasize that, as matters stand, our results have not lead to a 
comprehensive treatment of issues like exact integrability or quantization in the full, two- 
polarization case. Rather, they provide new windows to tackle these issues. In particular, 
the techniques introduced in Sees. Ill and V should enable one to analyze the reduced 
Einstein system along lines that are quite different from the traditional ones. 

Note added: After this work was completed and posted on the LANL archives (gr/qc 
9712053), we were informed of a paper by Romano and Torre [Q which contains a result 
which is equivalent to that presented in Sec. III.B. Their main interest is in the 'issue of 
time' and they work in a Hamiltonian framework with a phase space which is enlarged to 
incorporate a 'clock degree of freedom'. As a side-remark, they point out that their Hamilto- 
nian description can be derived from a symmetry reduced three-dimensional harmonic map 
in the parametrized field theory formalism. 
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